Geometry of phase separation 
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We study the domain geometry during spinodal decomposition of a 50:50 binary mixture in 
two dimensions. Extending arguments developed to treat non-conserved coarsening, we obtain 
approximate analytic results for the distribution of domain areas and perimeters during the dynam- 
ics. The main approximation is to regard the interfaces separating domains as moving independently. 
While this is true in the non-conserved case, it is not in the conserved one. Our results can therefore 
be considered as a 'first order' approximation for the distributions. In contrast to the celebrated 
Lifshitz-Slyozov- Wagner distribution of structures of the minority phase in the limit of very small 
concentration, the distribution of domain areas in the 50:50 case does not have a cut-off. Large struc- 
tures (areas or perimeters) retain the morphology of a percolative or critical initial condition, for 
quenches from high temperatures or the critical point respectively. The corresponding distributions 
are described by a cA~ T tail, where c and r are exactly known. With increasing time, small struc- 
tures tend to have a spherical shape with a smooth surface before evaporating by diffusion. In this 
regime the number density of domains with area A scales as A 1 ^ 2 , as in the Lifshitz-Slyozov- Wagner 
theory. The threshold between the small and large regimes is determined by the characteristic area, 
A ~ [\(T)t] 2 / 3 . Finally, we study the relation between perimeters and areas and the distribution of 
boundary lengths, finding results that are consistent with the ones summarized above. We test our 
predictions with Monte Carlo simulations of the 2d Ising Model. 



I. INTRODUCTION 

Phase separation is the process whereby a binary mix- 
ture of components A and B, initially in a homogeneous 
phase, demix leading to the coexistence of two phases, 
one rich in A and the other in B. The system, initially 
in an unstable spatially uniform state, performs a coars- 
ening process to approach its thermodynamically stable 
phase-separated state 

In the case of fluids, hydrodynamic effects may be im- 
portant to the demixing process Qj; these interactions 
are hard to treat analytically and the results of numeri- 
cal simulations are sometimes controversial. Some stud- 
ies even claim that dynamical scaling is broken by hy- 
drodynamic transport In the following we focus on 
the phase separation process in systems without hydro- 
dynamics. A typical realization is phase separation in 
binary alloys, high-viscosity fluids and polymer blends. 

The time-dependent order parameter characterising 
the phase separation phenomenon is a continuous scalar 
field, <p(f,t), that represents the local difference in con- 
centration of the two phases normalized by the sum of 
the averaged concentrations. Its evolution is described 
by a phenomenological Langevin-like equation: 
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H[(j>] is a Ginzburg-Landau type 'free-energy' with an 
clastic term and a double-well potential. The stochastic 
field r)(r,t) represents a conserved thermal noise, while 
the mobility M(<j>) is in general a function of the field <j>, 
as indicated. When M(dj) = 1, Eq. JT]) is known as the 
Cahn-Hilliard equation [4[ , or model B in the Hohcnberg- 
Halperin classification of critical dynamics [f|. Several 



discrete models to study the phase-separation process in 
binary alloys have also been proposed and studied in the 
literature. These are lattice gases that are themselves 
mapped onto Ising models with locally conserved order 
parameter 0, 0] ■ The dynamics follow stochastic rules for 
the interchange of nearest neighbour A and B molecules 
or, in the spin language, the reversal of a pair of neigh- 
bouring antiparallcl spins. 

Two microscopic processes contribute to phase order- 
ing dynamics with locally conserved order parameter, 
namely bulk and surface diffusion. In bulk diffusion a 
molecule separates (evaporation) from the surface of a do- 
main, diffuses within the neighbouring domain of the op- 
posite phase, and finally attaches to its original domain or 
another one. In the context of an Ising model simulation, 
there is an activation energy for this process and one can 
check that the dominant growth mechanism is the trans- 
port of material through the bulk from domain bound- 
aries with large curvature to domain boundaries with 
small curvature. In surface diffusion molecules 'walk' on 
the interface at no energy cost. This mechanism leads to 
the motion of whole domains in the sample and thus the 
possibility of merging two domains together when they 
collide. In the usual Kawasaki J6| spin-exchange dynam- 
ics or in model B's evolution 8], the early dynamics is 
surface diffusion driven and at later time the dominant 
process becomes bulk diffusion. There exist, however, 
modified versions where either bulk or surface diffusion 
are suppressed d, d, [l(| ■ 

For deep temperature quenches, M(4>) — ► 0, bulk dif- 
fusion is effectively eliminated, and domain growth pro- 
ceeds by surface motion. In quenches to moderate sub- 
critical temperatures, on the other hand, the mobility 
does not play an important role, M (</>) ~ const., and the 
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domain growth is bulk-driven. 

Let us now review in some detail the main features of 
coarsening with locally conserved order parameter. Using 
the Cahn-Hilliard equation, it can be easily shown that 
the radius, R(t), of a single spherical domain of negative 
phase (cf) = —1) in an infinite sea of positive phase (0 = 
+1), evolves from time t = to time t as 

R 3 (t) = R 3 (0)-~at. (2) 

with a a parameter that quantifies the surface tension. A 
domain with initial radius Ri thus evaporates in a time 
t ~ R 3 in contrast to the non-conserved order parameter 
dynamics in which the area within a boundary simply 
shrinks under the curvature force in a time t ~ -R? [lH. 

Lifshitz-Slyozov [l2j and Wagner [l3| studied, for a 
three-dimensional system, the growth and shrinkage of 
domains of one phase embedded in one large domain 
of the other phase in the limit of small minority phase 
concentration, c — > 0. In their celebrated papers they 
realized that domain growth at the late stages is lim- 
ited by matter diffusion through the majority domain. 
In this case the evolution of a domain of the minority 
phase with radius i?, immersed in a sea of the majority 
phase that is 'supersaturated' with the dissolved minority 
species [l], [l2| can follow two paths: the domain evapo- 
rates by diffusion if Ri < R Cl or it grows by absorbing 
material from the majority phase if Ri > R c , where R c 
is a time-dependent 'critical radius'. This critical radius 
turns out to be the only characteristic length-scale in the 
system, R(i), and serves to scale all correlation functions 
according to the dynamic scaling hypothesis. It grows as 
R(t) = R c (t) ~ t 1 / 3 . 

Lifshitz-Slyozov [13] also derived an expression for the 
density of droplets of the minority phase with linear size 
R in d = 3. Three important properties of the Lifshitz- 
Slyozov distribution are: 

• The distribution of droplet radii has an upper cut- 
off, -Rmax(i), where -R max (£) ~ t 1 ^ 3 is a constant, 
equal to 3/2, times the critical radius R c (t). 

• The decay close to the cut-off is exponential. 

• The density of small objects, R ~ 0, satisfies scal- 
ing and behaves as n(R,t) ~ R c (t)~ 4 (R/ R c (t)) 2 , 
where n(R, t)dR is the number of droplets per unit 
volume with radius in the interval (R, R + dR) . 

Later simulations established that the sca ling functions 
depend on the minority concentration fLU Ha ]. 

The Lifshitz-Slyozov calculation can easily be extended 
to any space dimension d > 2 [la) . However, the limit 
d — > 2 is singular, and does not commute with the limit 
c — * 0. Rogers and Desai [TtJ showed, however, that the 
usual scaling forms apply in d = 2, with R(i) ~ t 1 / 3 , for 
large t at small non-zero volume fraction c. 

More recently, Huse used scaling and energetic argu- 
ments to generalize the Lifshitz-Slyozov growth law and 



argued that it should also apply to critical quenches with 
equal volume fractions of the two phases [18| ■ Numerical 
simulations [H, 0], m agreement also with the earlier 
ones in Q , suggest that the typical domain radius scales 
in time as R(i) ~ t 1 / 3 for any value of c, even in the 50:50 
case. The scaling function for the distribution of domain 
areas has not been analyzed in this case. 

If we take also into account the competition between 
bulk and surface diffusion the growth law is modified at 
early times. The former process is the one responsible 
for the scaling of global observables with a typical do- 
main length R(i) ~ t 1 ^ 3 while the latter yields a slower 
time-dependence, R(i) ~ t 1 / 4 , that is important only 
at relatively short times after the quench [2(| HJ, l22j |. 
The temperature-dependent crossover can be seen, for ex- 
ample, in numerical simulations with Kawasaki dynam- 
ics [8(. The crossover time diverges when T — > 0. This 
observation has been used to develop accelerated algo- 
rithms to simulate discrete models in which only bulk 
diffusion processes arc considered, which should describe 
phase-separation correctly at late times after the quench. 

Phase-separation in the Kawasaki spin-exchange dy- 
namics is equivalent to a Cahn-Hilliard equation with or- 
der parameter dependent mobility. In 10], a model with 
M((f>) = 1 — a<f> 2 was studied. The time dependent struc- 
ture factor exhibits dynamical scaling, and the scaling 
function is numerically indistinguishable from the Cahn- 
Hilliard one, consistent with what was expected from nu- 
merical studies with Kawasaki dynamics. 

In this paper we study the morphology of domain and 
perimeter structures in the spinodal decomposition of a 
two dimensional system with equal concentrations of the 
two phases. In particular, we analyse the distributions 
of the domain areas and their associated perimeters, and 
the relation between areas and perimeters during the evo- 
lution. We consider bulk-and-surface diffusion and just- 
bulk diffusion processes. 

Extending a formalism previously developed for the 
study of domain growth in the non-conserved case [23l . 
HH Ell 1 we propose an analytic form for the domain size 
distribution function in its full-range of variation, and 
we test it with simulations on the two-dimensional Ising 
model (2dIM). Our analytic prediction for the distribu- 
tion of small areas is the result of one hypothesis: that 
interfaces move independently. This assumption is valid 
for domains of any size in curvature driven dynamics: the 
fission of a big domain into two smaller ones or the coa- 
lescence of two domains to form a bigger one are forbid- 
den in the continuous Allen-Calm description [ll[, and 
are not important in the heat-bath dynamics of the Ising 
model with non-conserved order parameter. In the con- 
served order parameter case, even with c — > 0, this as- 
sumption does not strictly hold and corrections must, in 
principle, be included. Indeed, already Lifshitz and Sly- 
ozov made an attempt to go beyond their simple model 
and account for coalescence when c — ► 0. Later, it be- 
came clear that in locally conserved order parameter dy- 
namics, the dominant effect not accounted for in the sim- 



3 



pie description that takes domain-boundaries as indepen- 
dent objects was interdomain correlations rather than co- 
alescence. For a discussion of the limit c — > considering 
interactions between droplets see (26[. 

Therefore, our analytical results are just a 'first order' 
approximation. Still, as we shall see in the main body 
of the paper, this approximation yields a very good de- 
scription of numerical data obtained with Monte Carlo 
simulations. The quantities on which we focus are do- 
main areas and 'hull enclosed areas', where the latter 
are the areas enclosed within the outer boundaries (or 
'hulls') of domains, and domain perimeters. The distri- 
bution functions for the hull enclosed and domain areas 
are nh t d(A, t), where rih,d{A, t)dA is the number of hulls 
(h) or domains (d) per unit area of the system with area 
in the range (A, A + dA) . 

The main properties of these distributions are: 

• The number density of domain and hull- 
enclosed areas satisfy scaling: nh,d(A, t) — 
t~ 4 ^ 3 fh,d(A/t 2 / 3 ). The argument of the scaling 
functions arises from the fact that the characteris- 
tic area of hulls and domains grows as t 2 / 3, . The 
prefactor i~ 4 / 3 follows from the fact that there is 
of order one domain per scale area. 

• The scaling functions, fh,d(x), do not have any cut- 
off and extend to infinite values of x falling-off as 
(2)cf li d£~ r , with x = A/ (\h,dt) 2 / 3 '• The constants 
Ch,d are the ones in the initial (or quasi-initial - see 
below) area distribution, rih,d{Ai) ~ (2)ch,d/AJ, 
for Ai — > co. The prefactor is known exactly, 
Ch = l/8ny/3 [27|. The factor 2 is present when 
quenching from high temperature, To > T c , and is 
absent when the initial condition is the critical Ising 
one, To — T c . The exponent t depends on whether 
we consider hull-enclosed (r = 2) or domain areas, 
and in the latter case on the initial condition, i.e. 
infinite or critical temperature. In both cases it is 
very close to 2. 

• After a quench from high temperature to a suf- 
ficiently low working temperature, T <C T c , the 
small- argument behaviour of the scaling function 
is f(x) oc y/x, in agreement with the Lifshitz- 
Slyozov- Wagner prediction for the small concentra- 
tion limit. At higher working temperature and for 
critical initial conditions the behaviour is modified 
in a way that we describe in the text. 

The paper is organised as follows. In Sect. HI] we de- 
scribe an approximate analytic derivation of the time- 
dependent hull enclosed and domain area distributions. 
These arguments do not rely on any scaling hypothesis 
but rather support its validity. We compare our a ppr oach 
to the celebrated Lifshitz-Slyozov- Wagner theory [12l [l3l] . 
In Sect. Hill wc show our numerical results for the statis- 
tics of areas in the 2dIM evolving with locally conserved 
dynamics. We use variants 9] of the Kawasaki algo- 
rithm Q that we briefly explain in this Section. The use 



of different algorithms allows us to switch surface diffu- 
sion on and off, and pinpoint the relative importance of 
these processes. Section [IV] is devoted to the analysis, 
both analytical and numerical, of the geometry of do- 
main walls during the dynamics and their relation to the 
corresponding areas. Finally, in the Conclusion we sum- 
marize our results and we discuss future studies along 
these lines. 



II. STATISTICS OF AREAS: ANALYTIC 
RESULTS 

In this Section we analyze the number density of hull- 
enclosed and domain areas. A domain is a region of 
connected aligned spins. Each domain has one external 
perimeter which is called the hull. The hull-enclosed area 
is the total area within this perimeter, i.e. the domain 
area plus the area of any internal sub-domains. 



A. Initial distribution 

We study the coarsening dynamics after a quench from 
To — * oo and To = T c . The first case is very similar 
to critical percolation, as we explained in [23| . Equilib- 
rium infinite temperature initial conditions (fraction of 
up spins = 1/2) are below the critical random percola- 
tion point (p c « 0.59) for a square lattice in d = 2. After 
a few MC steps, however, the system reaches the critical 
percolation condition, e.g. the expected A~ 2 tail is ob- 
served in the distribution of hull enclosed areas. We have 
checked that this is so from the analysis of several corre- 
lation functions as well as the distribution of structures. 
The same happens for non-conserved order parameter dy- 
namics [23|]. This fact justifies the use of the Cardy-Ziff 
exact result for the distribution of hull-enclosed areas at 
critical percolation [27j as our effective initial condition 
from To — > 00. From a more general perspective, one can 
argue that if the system is coarse-grained to the scale 
R(t), and R(t) is large compared to the lattice spacing, 
a continuum description becomes appropriate, for which 
the percolation threshold is precisely 1/2. At the scale 
R(i), the system is disordered and has the character of 
the percolation model at threshold. 

In Section II in [23[ we recalled the equilibrium distri- 
bution of hull-enclosed and domain areas, domain walls 
and their geometrical relation to their associated areas 
at critical percolation and critical Ising initial conditions. 
In Section III we derived some generic results that follow 
from the use of sum rules and the scaling hypothesis. We 
do not repeat the description of these properties here but 
we refer the reader to [23j for further details. 

For future reference we simply list here the initial dis- 
tributions, rih.d{Ai, ti), of hull-enclosed and domain areas 
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where the exponent r describes the critical Ising state 
and t' the critical percolation state. 

The expression for rih(Ai,ti) is asymptotically exact 
for areas Ai large compared to the microscopic area, Ao ps 
a 2 , where a is the lattice spacing. The factor 2 is present 
for critical percolation but absent for the critical Ising 
initial condition [27[ . In the former case, the "initial" time 
ti is taken to be a few MC steps, such that the system 
has effectively reached the percolative critical state. The 
constant Ch is known analytically, Ch = 1 /8nV3 [13] ■ 

The expressions for rid(Ai,ti) require a little more dis- 
cussion. The key point is that the constant Ch is very 
small (ch — 0.023), so that "domains- within-domains" are 
rare. As a consequence, the domain area distributions are 
very similar to the corresponding distributions of hull en- 
closed areas. However, as discussed at length in (23[, the 
domain area distributions must fall off faster than A~ 2 
in order that the total domain area, per unit area of the 
system, be finite (and equal to unity). The smallness of 
Ch suggests that they will fall off only slightly faster than 
A~ 2 , and indeed the exponents r and r' in Eq. take 
the numerical values r ~ 2.027 and r' ~ 2.055. For- 
mally, if Ch is regarded as a small parameter, one finds 
r = 2 + O(ch), t' = 2 + 0(c h ) and c d = c h + 0(c\) 
This line of reasoning also accounts for the factor 2 that 
(again to leading order in Ch) appears in the third of 
Eqs. ([3]) but not the second. 



B. Characteristic domain length 

It is by now well-established [H, EH that the spatial 
equal-time correlation function in demixing systems is 
correctly described by the dynamic scaling hypothesis 
with a characteristic length R(i) ~ t 1 / 3 in agreement 
with the Lifshitz-Slyozov- Wagner prediction jl2l fl3T] and 
the extension beyond the small concentration limit de- 
rived by Huse [18| . The effect of temperature fluctuations 
is expected to be described by a T-dependent prefactor, 
R(t, T) = [\(T)t]^ 3 . In the following we do not write 
the T dependence explicitly. 



C. Large structures 

It is natural to assume that at time t, as for non- 
conserved order parameter dynamics, large structures, 
characterized by a long linear dimension R R(i) , have 



not changed much with respect to the initial condition. 
We shall support this claim with the numerical results. 
Thus, for sufficiently large hull-enclosed areas such that 
the time-dependence can be neglected we expect 



n h (A,t) 



(2)c h 
A 2 



A^t 



2/3 



(4) 



Similarly, for large domains, the area dependence of 
their distribution follows that of the initial condition, 
Eq. Q. 



D. Small structures 

Small structures, such that R <C R(i), are mostly em- 
bedded in very large domains. To a first approximation 
we shall assume that they are independent. Moreover, 
they are not expected to have holes of the opposite phase 
within, implying the equivalence between hull-enclosed 
and domain areas at these scales. Indeed, any smaller 
structure placed within must have evaporated by time t. 
We then propose that the number density of small hull- 
enclosed or domain areas at time t can be written as a 
function of the initial distribution, 



n(A,t) 



dAi 5(A- A(t,Ai))n(Ai,ti) , (5) 



with Ai the initial area and n(Ai, ti) their number density 
at the initial time U. A(t,Ai) is the area of a domain 
at time t that had area Ai at time ti. In writing this 
equation we have implicitly assumed that an area cannot 
split into two, and that two such areas cannot coalesce, 
which is not strictly true in conserved order parameter 
dynamics. 

Note that for sufficiently large areas so that the time- 
dependence is not important and A ~ Ai Eq. ([5]) immedi- 
ately yields n(A, t) = n(Ai,ti) and Eq. ([J| is recovered. It 
has to be stressed, however, that Eq. ([5]) does not strictly 
apply in this case. 

Assuming that the small areas are circular 
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see Eq. $2$ , and after a straightforward calculation using 
Eq. ^ for the initial distribution one finds 
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for hull-enclosed areas. This prediction has the expected 
scaling form nh(A, t) — i~ 4 / 3 f(A/t 2 / 3 ) corresponding to 
a system with characteristic area A(t) ~ t 2 ! 3 or charac- 
teristic length scale R(t) ~ i 1 / 3 . At very small areas, 



5 



A <C (\ht) 2 / 3 , where our approximations are better jus- 
tified, one has 



(X h t) 4 ' 3 n h (A,t) n (2)c H 



A 



1/2 



(8) 



As expected, taking the limit A/(X h t) 2/3 > 1 in Eq. Q 
one recovers Eq. (|3]) . Although this limit goes beyond the 
limit of validity of Eq. |0 we shall propose that Eq. (|8|) 
actually holds, at least approximately, for all values of 

In [231 ] we studied non-conserved order parameter dy- 
namics and we derived the number density of domain ar- 
eas from the one of hull-enclosed areas. The key fact in 
this case was that we could treat the distribution of hull- 
enclosed areas exactly and then use a small Ch ~ 0.023 
expansion to get the statistical properties of domains. In 
the case of phase separation our results for hull-enclosed 
areas are already approximate. Still, the relation between 
hull-enclosed area distribution and domain area distribu- 
tion obtained in [23| should remain approximately true, 
as a first order expansion in small ct , since small domains 
are not expected to have structures within. In conclusion, 
for large areas and long times such that a regularizing 
microscopic area Aq = X d to can be neglected, we expect 
the same functional form as the one given in Eq. CZJ with 
c d = c h + 0(cl), X d = X h [l + 0(c h )] and the power 5/3 
in the denominator replaced by (2t+ l)/3 or (2r' + l)/3: 



(X d t) 4 / 3 n d (A,t) 



= g(x) 

with x = A/(X d t) 2 / 3 , and 
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(X d t) 4 / 3 n d (A,t) ~ (2)c d 
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(10) 



for A -C (X d t) 2 / 3 . This expression can be compared to 
Eq. (49) in [23| valid for non-conserved order parameter 
dynamics. 

Let us emphasize again the main approximation of our 
analytical approach: we are considering each domain area 
as an independent entity. This is strictly true for the 
non-conserved order parameter [23l . |24| but is an approx- 
imation in the conserved order parameter problem. 



dynamical universality class, the scaling functions being 
the same, then Eq. |(7J) can be generalized in the form 
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R 2 (i) 



and similarly for the domain area distribution. The time- 
dependence in R(t) should include all regimes (e.g. t 1 / 4 
and i 1 / 3 in the clean case with surface and bulk diffusion) 
and can be extracted from the dynamic scaling analysis 
of the correlation functions. We shall check numerically 
the super-universality hypothesis. 



III. STATISTICS OF AREAS: NUMERICAL 
TESTS 

To test our analytic results we carried out numerical 
simulations on the 2d square-lattice Ising model (2dIM) 
with periodic boundary conditions. 



H 



(12) 



where cr, = ±1. We will start considering the pure 
model, Jij — J > Vij, and then the random bond model 
in which the Jij are random variables uniformly dis- 
tributed over the interval [1/2, 3/2]. We used several ver- 
sions of conserved order parameter dynamics that switch 
on and off surface diffusion. These are Kawasaki dynam- 
ics at finite temperature including both surface and bulk 
diffusion, and accelerated bulk diffusion in which surface 
diffusion is totally suppressed. Bulk diffusion needs to 
overcome energy barriers; thus this variant runs at fi- 
nite temperature only. In all cases we implemented the 
continuous time method and the algorithms become re- 
jection free. A detailed description of these algorithms 
appeared in [9j. Domain areas are identified with the 
Hoshen-Kopelman algorithm [30l ] . 

All data have been obtained using systems with size 
L 2 = 10 3 x 10 3 and 10 3 runs using independent initial 
conditions. We ran at different temperatures specified 
below. We considered two types of initial conditions, 
equilibrium at infinite temperature, To — > oo, and equi- 
librium at the critical Ising temperature, To — T c . All 
results are presented in equivalent Monte Carlo (MC) 
times. An example, for To — > oo, is presented in figure [TJ 
with snapshots taken at different times. 



E. Super universality 

All the results above are valid for the bulk diffusion 
driven case. What happens if we consider the case in 
which bulk and surface diffusion are in competition, or 
whether we include quenched disorder in the couplings? 
If we suppose that all these systems belong to the same 



A. Domain areas 

In this Section we show data obtained with the accel- 
erated bulk diffusion algorithm but we checked that simi- 
lar results are obtained with finite-temperature Kawasaki 
dynamics. 
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FIG. 1: Snapshots of the 2d Ising model evolving with locally 
conserved order parameter at t = 1000 and t = 16000 equiva- 
lent MC steps using the accelerated bulk algorithm explained 
in Sect. Hill The initial condition is a random configuration 
of ±1 spins taken with probability one half, that is to say an 
infinite temperature state. Evolution occurs at T = 1.0. 



1. Infinite temperature initial condition, low T evolution 

In Fig. [3 (top) we show the time-dependent domain 
area distribution in double logarithmic scale, at three dif- 
ferent times, following a quench from To — ► oo to T = 1.0. 
The working temperature is very low compared to the 
critical value, T c = 2.269. 

The figure shows a strong time dependence at small ar- 
eas and a very weak one in the tail, which is clearly very 
close to a power law. The curves at small areas move 
downwards and the breaking point from the asymptotic 
power law decay moves towards larger values of A for 
increasing t. We include the spanning domain in the 
statistics: the bumps on the tail of the distribution is a 
finite size effect visible only when the number of domain 
areas has already decreased by several orders of magni- 
tude. In the tail of the probability distribution function 
(pdf) the numerical error is smaller than the size of the 
data points. The discussion of finite size effects in [23[ 
also applies to this case. 

We test the analytic prediction: the very good agree- 
ment between the analytical theory and the data is quite 
impressive. In Fig. [2] (bottom) we scale the data by 
plotting (X d t) 4 / 3 n d (A, t) against A/(\ d t) 2 / 3 with X d = 
0.0083. For A larger than the 'typical' value (\ d t) 2/3 
the time and X d dependence become less and less im- 
portant. We fit the parameter X d (T) by analysing the 
behaviour at small areas, A 3 / 2 < \ d {T)t, and we find, 
that \ d (T = 1.0) = 0.0083 yields the best collapse of 
data. We use the value c d = 0.025 (23[. The full line is 
our prediction Eq. (|9"|). 

In Fig. [3] we present the domain area distribution for 
the evolution using the Kawasaki algorithm both for the 
clean and the disordered system. The agreement with the 
analytical prediction is as good as with the bulk-diffusion 
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FIG. 2: (Colour online.) Top panel: Number density of do- 
mains areas per unit area for the 2dIM evolving at T = 1.0 
after a quench from To — > oo using the accelerated algo- 
rithm. The dotted line represent the analytical prediction. 
Bottom panel: rescaled data using the typical domain area 
time-dependence A ~ t 2 ^ 3 . The dotted line is the theoretical 
prediction in Eq. © with r replaced by r', appropriate to a 
disordered initial state. 



algorithm, suggesting the validity of super universality 
between both dynamics. 

We extract the growing length R(i) from the analysis of 
the spatial correlation function [23| and we find very good 
agreement between the numerical data and the scaling 
function suggesting that super-universality with respect 
to the inclusion of disorder in the interactions also holds. 



2. Critical temperature initial condition, low T evolution 

We now use a critical Ising initial condition, To = T Cl 
and show the results in Fig. |4j In this case we have to 
further distinguish large and small domain sizes. We find 
that the distribution at large areas, A ^> (\ d t) 2 / 3 is well 
described by the initial condition form. At the crossover 
A ~ {\ d t) 2 / 3 the data points leave the asymptotic power- 
law with parameter c d to approach the one with 2c d . For 
small areas, A < (X d t) 2 ^ 3 we find that Eq. © is not 
satisfied, see Fig. 2] Although few data points fall in this 
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FIG. 3: (Colour online.) Simulations with the Kawasaki al- 
gorithm, without and with disorder. The solid lines repre- 
sent the analytic prediction with c<j = 0.025, r' = 2.055 and 
Ad = 0.3. For the disordered case, the growth law R(t) is ex- 
tracted from the spatial correlation function. Even if there is 
deviations from the theoretical curve, the data seems to show 
sup er- universality . 



regime the departure from y/x is clear. 



by T and all data are well-described by the initial con- 
dition form. Instead, the small scale behaviour depends 
strongly on temperature fluctuations. The anomalous 
up-rising part of the distribution at small areas is time- 
independent, suggesting that it can be associated with 
equilibrium fluctuations of the domain-walls. As shown 
in Fig. O it is possible to extract the interface thermal 
fluctuations as was done in ref. [23} with the bulk equi- 
librium droplets (see Figs. 20 and 21 of that paper). 
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FIG. 4: (Colour online.) Number density of domain areas per 
unit area for the 2dIM evolving at T = 1.0 after a quench 
from To = T c using the accelerated bulk algorithm. The lines 
represent Eq. ((9J with constant Cd and 2cd. 

We believe that the reason why the analytic prediction 
for small areas fails in this case is that the assumption 
of independence of small domains is not justified for a 
critical Ising initial condition. 

3. Dependence on the working temperature 

Up to now we showed results obtained using a rather 
low working temperature. We now study whether and 
how our results are modified when using higher values of 
T. Figure[5]shows the numerical data and proves that the 
scaling is well-satisfied at all times and for all areas. The 
large scale behaviour of the distribution is not modified 



FIG. 5: (Colour online.) Number density of domains areas per 
unit area for the 2dIM evolving at T = 1.5 after a quench from 
To — » oo using the accelerated bulk algorithm. The prefactor 
Xd(T) in the growing-length is chosen to be Ad(T) = 0.0060. 
Compare to the results shown in Fig. [2] obtained for a lower 
working temperature. 



IV. STATISTICS OF PERIMETERS AND 
FRACTAL PROPERTIES 

The analytic argument described in Sect. |TT] can be ex- 
tended to study the distribution of domain wall lengths or 
perimeters. The domain perimeter is the total length of 
the interface between the chosen domain and the neigh- 
bouring ones - including the hull and internal borders. 
In this Section we present the analytic prediction for this 
function together with numerical results that confirm it. 
We concentrate on To — > oo and low-working temper- 
ature. In the simulations we define the length of the 
boundary as the number of broken bonds. 

A detailed description of the relation between the do- 
main areas and their boundaries in the critical Ising and 
critical percolation conditions as well as the number den- 
sities of domain perimeters was given in [23|. Here we 
focus on their evolution. 

After a quench from To — > oo, the domain areas, A, 
and their corresponding perimeters, p, obey the scaling 
relations (see Fig. [6]) 
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The relation between areas and perimeters exhibits two 
distinct regimes with a quite sharp crossover between 
them. During the coarsening process a characteristic 
scale R(*) ~ (Ad*) 1 / 3 develops such that domains with 
area A > R 2 (*) have the same exponent as in the ini- 
tial condition (structures that are highly ramified with 
a f< ~ 1) and domains with A < R 2 (t) become regu- 
lar (a' ~ 2). Interestingly, the small structures in the 
non-conserved order parameter dynamics are not com- 
pletely circular, as demonstrated by the fact that their 
a'< ~ 1.8, see Fig. [6] 
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for small areas, Aj (Ad*) 2 / 3 < 10, and the same expression 
with r]' d and a' d replaced by rj' d and a' d for large areas 
A/(\ d t) 2 / 3 > 10. Note that these expressions satisfy scal- 
ing. The scaling function, /<(a;), with x — p/(X d t) 1 ^ 3 , 
reaches a maximum at 
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and then falls-off to zero as another power-law. There 
is then a maximum at a finite and positive value of p as 
long as a'/*' > 1, that is to say, in the regime of not too 
large areas. The numeric evaluation of the right-hand- 
side yields x max = p m ax/{^dt) 1/3 ~ 3 which is in the 
range of validity of the scaling function /< . The time- 
dependent perimeter number density for long perimeters 
falls-off as a power law />(x) - a; Q '^ 1 - T ')- 1 . Although 
the function /> also has a maximum, this one falls out 
of its range of validity. The power law describing the tail 
of the number density of long perimeters is the same as 
the one characterising the initial distribution. 

In Fig. [71 top and bottom, we display the time- 
dependent perimeter number densities for a system evolv- 
ing at T = 1.0 after a quench from To — > oo. Notice that 
the perimeter length definition that we use on the lattice 
can only take even values and thus when constructing the 
histogram we have to take into account the extra factor 
of 2 in the binning. The data are in remarkably good 
agreement with the analytic prediction; the lines repre- 
sent the theoretical functional forms for long and short 
lengths, and describe very well the two limiting wings of 
the number density. The maximum is located at a value 
that is in agreement with Eq. (fTT)) . 



FIG. 6: (Colour online.) The time-dependent relation be- 
tween the area and the domain boundary evolving at T = 1.0 
after a quench from To — > oo using conserved and non- 
conserved order parameter dynamics. The lines are fits to the 
data points. For the large structures the fit yields ct > ~ 1 
in both cases while for the small structures it yields a ~ 2 
for conserved order parameter dynamics (circular domains) 
and a' y ~ 1.8 for non-conserved order parameter dynamics. 
Here z = 2 and 3 for non-conserved and conserved dynamics 
respectively. 

In analogy with the derivation in Sect. [II] for the 
time-dependent number density of domain areas, the 
time-dependent number density of domain- wall lengths, 
n d (p,t), is given by 



(A d *) n d (p, t) 



. (Ad*) 3 



l + (^) 3/2 



(Ad*) 1 / 5 



(16) 



V. CONCLUSIONS 

In this paper we studied the statistics and geometry 
of hull-enclosed and domain areas and interfaces during 
spinodal decomposition in two dimensions. 

The analytical part of our work is an extension of what 
we presented in [23J for the non-conserved order parame- 
ter case. The numerical part of it deals with Monte Carlo 
simulations of the 2rfIM with locally conserved magneti- 
zation. Our main results are: 

(i) We derived the scaling functions of the number den- 
sity of domain areas and perimeters with an approximate 
analytic argument. The expression that we obtained has 
two distinct limiting regimes. For areas that are much 
smaller than the characteristic area, R 2 (*), the Lifshitz- 
Slyozov- Wagner behaviour is recovered after a quench 
from To ~~ * 00 an d evolving at sufficiently low T. These 
structures are compact with smooth boundaries, close to 
circular, since the area-perimeter relation is A ~ p 2 . 

(ii) At higher T the small area behaviour departs from 
the Lifshitz-Slyozov- Wagner prediction. As for non- 
conserved order parameter dynamics, once we subtracted 
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FIG. 7: (Colour online.) The time-dependent number density 
of perimeters evolving at T = 1.0 from an initial condition at 
To — > oo. Top panel: raw data; note that the time-dependence 
is visible in the whole range of values of p (while in the area 
number densities the large area tails were very weakly de- 
pendent on time). Bottom panel: scaled data and analytic 
predictions for the small and large area regimes. 



the contribution from thermal domains within the grow- 
ing structures, the universal prediction is recovered. 

(iii) For critical Ising initial conditions large structures 
keep the geometry they had initially, but small struc- 
tures fail to follow our prediction. We conjecture that 
the reason is that our starting assumption, independence 
of domain wall motion for small domains, is not valid due 
to strong correlations in this case. 

(iv) The geometrical properties and distribution of the 
time-dependent areas that are larger than R 2 (t) are the 
ones of critical continuous percolation (for all initial con- 
ditions equilibrated at To > T c ) and critical Ising (for 
To = T c ). The long interfaces retain the fractal geometry 
imposed by the equilibrium initial condition. 

These results complement the ones that we presented 
in [H, [25[ . It would be interesting to experimentally 
check these results, as was done for non conserved dy- 
namics in [3TI ] . It is important to summarize the picture 
that emerges from all these studies. 

Large structures, with area larger than the typical one 
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FIG. 8: (Colour online.) Comparison between the number 
density of domain areas in the clean and disordered 2d Ising 
model with non-conserved dynamics and conserved dynamics 
in a zoom over the small area regime. The data points are the 
result of the numerical simulations while the straight lines are 
the analytic predictions. 



A > R 2 , keep the geometry they had initially in all re- 
spects. The scaling functions are thus independent of the 
dynamics - conserved or non-conserved - and the pres- 
ence or not of weak disorder. 

Small structures do not keep the initial geometry but 
become more circular in all cases. For a clean system 
with conserved order parameter, domains are indeed cir- 
cular while for non-conserved order parameter they be- 
come close to circular. The scaling function of their dis- 
tribution depends on the type of dynamics considered but 
it does not depend on the presence of weak disorder once 
scaling by the pertinent growth law has been taken into 
account. In Fig. [5] we illustrate this statement by com- 
paring the number densities of domain areas, on the one 
hand, in the random-bond Ising model (RBIM) with non- 
conserved order parameter (red datapoints) and the ana- 
lytic result for non-conserved order parameter dynamics 
and, on the other hand, in the clean and disordered con- 
served order parameter dynamics. 

Once the super-scaling of the number density with re- 
spect to weak quenched randomness is established numer- 
ically we can infer what it implies for the time-evolution 
of the individual small structures under the assumption 
of independence discussed in the text. Indeed, starting 
from Eq. © one can derive Eq. (fTTj) if one assumes that 
each time-dependent area is linked to the initial value by 



A z l 2 (t) « A, z/2 -R(t) 



(18) 



where z = 2 for systems in the universality class of model 
A (e.g. the non-conserved order parameter dynamics of 
the RBIM), and z = 3 for systems in the universality 
class of model B (e.g. the Kawasaki dynamics). Note 
that inside each universality class, the growth laws R(£) 
can be very different. 
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